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Abstract 

Charge and spin response in the planar t — J model at finite temperatures 
are investigated numerically, in the regime corresponding to cuprates at in- 
termediate doping. We show that the local spin correlation function is T- 
independent, leading to x"(u; < J) oc tanh(w/2T). The current correlation 
function C(u) appears to be T- as well as w-independent, and the optical con- 
ductivity a oc (1 — e~ w / r )/tj. Such anomalous response functions are brought 

in connection with the large entropy persisting at T < J. 
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Normal-state properties of cuprates offer a clear message, that metals with strongly 
correlated itinerant electrons are not yet properly understood. Although the main goal 
remains to relate the correlation effects with the onset of superconductivity at high T c , 
anomalous static and dynamical quantitites at T > T c alone represent a challenging subject. 
From recent experiments testing normal state properties such as the dc resistivity |I| , optical 
conductivity [@|§, neutron scattering [|]||, and the NMR and NQR relaxation ||, a unifying 
picture seems to emerge Q. These properties mainly depend on the extent of doping, and 
the regimes have been classified into the underdoped, the 'optimal' (intermediate-doping) 
and the overdoped, respectively. 

We investigate here the intermediate-doping (optimal) regime. In experiments the latter 
is usually defined by highest T c , whereas we refer only to the characteristic normal-state 
properties (e.g. the linear p(T) law), as summarized within the marginal Fermi liquid (MFL) 
hypothesis 0. Consistent with experiments, the low-frequency dynamics in cr(u) has been 
modeled within the Drude form by an anomalous effective relaxation rate r _1 oc to + r]T, 
and the spin susceptibility with ~x"{oj < T) oc uo/T. Alternative phenomenological theories 
|8||5[] have been also proposed, nevertheless the origin of these phenomena is not yet settled. 

We discuss finite-T (normal-state) properties of doped antiferromagnets (AFM) within 



the t — J model [10 



H = -t {c) s c is + H.c.) + J XX<Si ■ Sj - -riiTij), (1) 

(ij)s (ij) 

where c\ s (c is ) are projected fermionic operators, prohibiting double occupancy of sites. The 
ground state of the t — J model and related properties have been intensively studied both 



by analytical [plj] and numerical methods [11], nevertheless some basic questions, e.g. con 



cerning possible pairing in the ground state, remain unsolved. 



T > properties have been approached by the high-temperature series expansion ||lz 
and recently by the present authors using a novel numerical method, based on the Lanczos 
diagonalization method combined with a random sampling [13]. The latter method has 



the advantage that it allows the study of dynamical response functions within the most 



challenging regime T,u < J, where one could hope for a universal behavior, as observed 



in the experiments on cuprates. Results concerning the charge [T^] and spin [HJ response 
in the intermediate-doping regime, obtained by this method, can be summarized as follows: 
(i) <j(uj) shows a non-Drude fall-off, consistent with the a»-dependent relaxation rate within 
the MFL concept |7[] and in cuprates (ii) qualitative as well as quantitative results 

for <j(uj) and p(T) agree reasonably well with the experimental ones, (iii) dynamical spin 
susceptibility x"(<f, u;) shows the coexistence of the high-frequency (u oc t) free-fermion- 
like contribution and the low-cu spin- fluctuation contribution, (iv) x"(q> u < T) shows a 
pronounced T dependence, consistent with the MFL form and (even quantitatively) with 
the NMR relaxation in cuprates. 

Above results indicate that the t — J model represents a promising framework for the 
study of the anomalous metallic state in cuprates. In this Letter we present the evidence for 
some important features, which can lead to a more consistent picture of the intermediate- 
doping regime: (i) local spin (temporal) correlations appear to be particularly universal 
and thus fundamental to the understanding of x(<L u ) response, (ii) conductivity <j{uj) is 
consistent with the fast-decaying T-independent current correlations, leading to a universal 
form for (t(u), representing an alternative to the MFL Ansatz, (iii) both facts seem to be 
in a close connection with the large entropy persisting down to low temperatures T <C J in 
the 'optimal' regime. 

Let us consider the dynamical spin response as given by the susceptibility x(<?> u ) an d 
the corresponding dynamical spin correlation function S(q, uS) 

X ''(g>) = (l- e -^)%>), 

roo 

S(q,cu) = Rej o dt e^(Si(t)S z _ ? ), (2) 

with (3 = 1/T (we use units with % — Ub — 1)- We first analyse the local spin correlation 
function Sl{ui) and its symmetric part S{ui) 
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S L ^) = -Y,S(q 



UJ 



I 5 } ? 

<7 



S(u) = S L (u) + S L (-u) = (1 + e-P»)S L (u). (3) 



It should be pointed out that Sl(uj) and the related susceptibility Xl(o>) are directly mea- 
sured (in cuprates) by neutron scattering @||. The NMR relaxation as well yields the 
information on Sl(u — > 0) |H§, provided that the AFM spin fluctuations q ~ Q = (n, n) 
are dominant. 

An important restriction for S(u) is the sum rule 

S(u)duj = n((Sn 2 } = j(l-c h ), (4) 
o 4 

where = Nh/N is the hole concentration. Since in any finite system S(q = 0,u) is ill- 
defined (due to conservation of total S z ), the q = term is omitted in Eq. (0), and the sum 
rule (f|) serves as a useful test. 

We perform the evaluation of S(u) via Eq. (^|) by calculating S(q, uj) using the finite-T 
diagonalization method for small systems, in this case for the t — J model on the square 
lattice with N = 16 — 20 sites. We fix J/t = 0.3 to remain in the regime of cuprates [JTOj . 



For the description of the method we refer to Ref. |13[], its application to the spin dynamics 
is given in Ref. ||15| . We stress again that the method yields macroscopic-like results for 
T > T*, while below T* finite size effects become pronounced. Within the intermediate 
regime 2/16 < Ch < 4/16 we find typically T* ~ 0.1 t. T* is larger both within the 
underdoped and the overdoped region (at fixed N). 

In Fig. 1 we display the S(u) for Ch = 1/20,3/16 and several T in the range 0.1 < 
T/t < 0.7. It is immediately evident that S(uj) at 'optimal' doping Ch = 3/16 is essentially 
T-independent in a wide T range, although one crosses the exchange-energy scale T ~ J. 
For the underdoped case Ch = 1/20 the behavior is analogous for higher T > T ~ 0.7 J 
(the same holds for the undoped AFM, and to some extent for Ch = 2/16), consistent with 
the quantum critical regime within the AFM Deviations at lower T < T (where the 



renormalized classical regime |I6] is expected in the AFM) could be an indication for the 
onset of a 'pseudogap', but we cannot exclude that these phenomena are finite-size artifacts 
since T > T*(c h ). 

To follow the doping dependence we present in Fig. 2 the variation with Ch at fixed 



T = 0.2 t < J. Here we plot the integrated intensity 

I s {u)= rS{u/)du', (5) 



since the latter does not require any smoothing. Again, we notice that for chosen T results 
are most reliable at intermediate doping, being poorer otherwise. The most striking message 
is that the initial slope of Is{w) and consequently Sl{oj — > 0) are nearly doping independent 
for < Ch < 0.25 (as well as T— independent at intermediate doping). This is consistent 
with the NMR (NQR) relaxation rates T{~ x as measured in La2- x Sr x Cu04 in the range 
x = 0-0.15 H[T5|. 

Only for the overdoped systems with > 0.25 the low-frequency behavior changes 
qualitatively, where the latter part is strongly supressed as expected in (more) normal Fermi 
liquids. Ig(u > J) is doping dependent even for < 0.25, consistent with the c^-dependence 
of the sum rule, Eq. (|). In addition, at the intermediate doping S(u) decreases smoothly 
(see Fig. 1) up to uj ~ 4t, this being the consequence of a free-fermion-like component [15| . 
On the other hand, in the underdoped regime the dynamics is restricted to uj < 3 J < t. 

How can one explain the universality of S(u) at intermediate cyP. First we note that up 
to Ch ~ 0.3 the dominant scale of spin fluctuations remains related to J. From the explicit 
expression in terms of the eigenstates of the system 

n,m 

x5(lu- E m + E n ), (6) 

one would conjecture (quite generally) a T-independence of response for uj ^> T. Although 
such plausibility arguments have been used previously, their validity clearly depends on the 
character and the density of low-lying many-body states (furtheron we present an evidence 
that the 'optimally' doped AFM has a particularly large density of low-lying states). To 
explain also the T-independence of S(u < T), we only need to recall the sum rule Eq. (^) and 
to assume that there is no characteristic scale uo c < T which could introduce an additional 
low-cu structure in S{uj). 



A natural scale for an AFM is the (gap) frequency uo c ~ c/£ where c is the spin wave 
velocity and £ the AFM correlation length. Here originates the essential difference between 
the undoped and the 'optimally' doped AFM. While for an AFM in the renormalized classical 
regime £ is exponentially large for T <^ J, and consequently u c < T in the doped case 
£ < 1/y/ch is determined predominantly by Ch, so £ is rather T- independent for T < J, 
excluding u c < T < J. 

We conclude the discussion of spin dynamics by consequences of the universality of S(ui). 
The local susceptibility is given by 

x£(«)=tanh(^)3( W ). (7) 

Since neutron scattering probes only u < J, one can simplify Eq. (0) further by S(u) ~ So. 
Such form, consistent with the MFL 0, has been recently used to describe experiments ||. 
Here, one should take into account that we are not able to establish within the t — J model 
the existence of the 'pseudogap' u g ~ 0.1 J ||, observed in cuprates at low T > T c . Further 
it follows from Eq. (0) that for T < J the relevant scale for x'l( u ) is ^ ~ 2T. The same 
should hold for the response at fixed q. So one can generalise Eq. (0) 

III-* \ Xq ii/ \ 27T In (^q-in) II, \ la \ 

X (q,w) -Xl(u) ~ ~ -5- -^XlH, (8) 

xl \q-Q\ +? 

where the cutoff q m ~ 7r. The scaling Eq. (|8|) is expected to hold only for |<f — Q| ^ 
£ _1 , where one should take into account that in the 'optimal' regime of the t — J model 
£ < 1 p^,I5|. Nevertheless, £ remains T-dependent (becoming even shorter at higher T), 



introducing additional T- variation in Eq. (|8|). Outside of the mentioned regime, in particular 
for q ~ 0, the response is more free fermion-like, i.e. \"{q,uj) is T-independent [ |I"5[ . 
Let us investigate in an analogous way the dynamical conductivity cf{uj) 

a{u) = C(uj), C(u) = Re / dt e^(j(t)j), (9) 



U) Jo 

where j is the current density (we put eo = 1). Unlike S(w), C{uj) does not obey a T- 
independent sum rule. Nevertheless, motivated by universal spin dynamics we reexamine our 
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results on a(u), obtained by the same finite- T numerical method on systems with N = 16—20 
fl4| . In Fig. 3 we present the corresponding integrated spectra Ic{w) = Jo C(u')du' for fixed 
doping in the 'optimal' regime Ch = 3/16, for various T < t. We establish several remarkable 
features: (i) for T < J spectra Ic{oJ) are essentially independent of T, at least for available 
T > T*, (ii) at the same time the slope of Ic(oJ < 2 t) is nearly constant, i.e. C{uo) ~ Co 
in a wide u;-range, Co being weakly J-dependent (tested for J ft = 0.2,0.6), (iii) even for 
higher T > J the differences, e.g. in the slope Co and in the sum rule Jc(oo), appear 
as less essential (note that for T ^> t we know exactly Jc(oo) = t 2 Ch(l — Ch) 13). We 
reproduce the same characteristic behavior also for Ch — 4/16 (confirming Co oc c^). On the 
other hand, Ch = 2/18 seems to represent a crossover to the underdoped regime, where in 
contrast (e.g. for = 1/20) we find a pronounced T- and independence of C(u>) at T < J. 
This qualitative change can be again attributed to a 'pseudogap' appearing at larger T in 
underdoped systems |]. 

Restricting our discussion to the intermediate doping, we note that C(lu) ~ C implies 
a nonanalytic behavior of a(u — > 0), starting with a finite slope at u = 0. This has been 



already evident in our previous results [13]. Moreover, with C(u) = Co we can claim a 
simple universal form for u < 2 t 

1 — 

a(u;) = Co . (10) 

UJ 

It is rather surprising that this form can be well fitted (for ui, T <C t) with a Drude-type 
form with an MFL effective relaxation rate r _1 = 2h\{uj + r]T) with specific A ~ 0.09 and 
7] ~ 2.7 [14j. The form Eq. (|10D for <t(co') trivially reproduces the remarkable linear law p oc T 



in cuprates [|IJ, as well as the non-Drude falloff at uo > T. A quantitative comparison of 
our <j(u) with experimental results [|2|,0] has been already established in Ref. ||14]| , including 



the MFL form for r. It is evident that the expression (|H]) is universal, containing only 
Cq as a parameter (the MFL Ansatz contains at least one more) and should be retested in 
more detail with experiments. More generally, for larger uo one should replace Co with an 
universal C(u), which could be however sensistive to a particular model, as could be also 



It should be mentioned that C(u) ~ Cq has been derived for a single hole conductivity 



within the retraceable path approximation [17]], with a restricted validity for T > t (or pos- 



sibly T > J). We find this behavior only for the intermediate doping, hence new arguments 
are needed. We can follow the analysis analogous to S(u>), Eq. expressing C(uS) in terms 
of eigenstates. As before, the T-independence of C(u > T) seems plausible. The fact that 
C(ui < 2t) ~ Cq requires however that the current relaxation is very fast, i.e. determined 
only by the incoherent hopping and the interhole collisions. The spin fluctuation scale J 
does not enter directly, i.e. even for T J the spin system only serves as a random bath 
for charge degrees of freedom (holons). This conclusion remains valid as far as there is no 
characteristic frequency u c < T (e.g. the 'pseudogap') in the system. 

Above arguments, both for Sl(u>) and C(u), require a large density (degeneracy) of 
the low lying many-body states, apparently being a crucial feature for the most challenging 
intermediate doping regime. To quantify this statement we calculate the entropy density s = 
S/N, using again the finite-T diagonalization method, being less space and time consuming 
for static quantities In Fig. 4 we present results obtained for N = 18 at various c^. Here 
Ch varies continuously, since we have to employ a grand canonical distribution. Our data 
agree with the high-T expansion results by B. Putikka (unpublished). The main lesson from 
Fig. 4 is that 'optimal' cases with ~ 0.1 — 0.3 are characterized by the largest entropy s at 
low T < J, e.g. s > 0.2/site for T ~ 0.2 J, e.g. being almost one half of s(T = oo) for the 
AFM. This implies a very large degeneracy of low-lying states, which could be attributed to 
the spin subsystem, frustrated by the hole motion. 

In conclusion, our results for the t — J model seem to indicate that the anomalous, but 
universal, dynamics in the intermediate-doping regime of correlated systems is a consequence 
of the extreme degeneracy and the collapse of low-lying quantum states (introduced by 
doping the magnetic insulator), which leads to a diffusive-like charge and spin response where 
T represents the only relevant energy (frequency) scale. On the other hand, our results do 
not exclude the possible onset of coherence (related to 'pseudogaps' and superconductivity) 
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at lower T <T*. 
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FIGURES 

FIG. 1. Local spin correlation function S(uj) for = 1/20, 3/16 and various T. 



FIG. 2. Integrated spectra is(u>) at fixed T = 0.2 t and various c^. 



FIG. 3. Integrated current correlation spectra Ic{w) for = 3/16 and various T. 



FIG. 4. Entropy density s (in units of fee) vs. T for various dopings c^. 
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